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Abstract We propose a Lemaitre-Tolman-Bondi system mimicking a two- 
body system to address the problem of the cosmological expansion versus 
local dynamics. This system is strongly bound but participates in the cosmic 
expansion and is exactly comoving with the cosmic substratum. 
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1 Introduction 

The problem of the influence of the cosmological expansion on local gravita- 
tionally bound systems was apparently first raised by McVittie in 1933 [I], 
studied by Einstein and Straus [21[3]i and then debated in dozens of papers 
stretching to our days (see the recent review by Carrera and Giulini [4 ). 

The problem is the following: for a local system such as a planetary sys- 
tem, a binary stellar system, or a galaxy in the expanding universe, does the 
cosmological expansion affect the local dynamics of this system and, if so, in 
which way and, numerically, to what extent? Many authors support the view 
that the cosmic expansion affects only systems larger than a certain spatial 
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scale and there is no effect below that scale. However, it has not been possible 
to assess what this scale is, or what determines it [HIS]. It seems that, assum- 
ing in principle the Friedmann-Lemaitre- Robertson- Walker (FLRW) metric to 
be valid down to small scales (possibly with some modifications that describe 
local inhomogeneities [IJinilZllH] ) ) all systems are subject to the effect of the cos- 
mic expansion, although this effect is numerically negligible for small systems 
and stronger for larger and larger objects, up to the scale of galaxy clusters 
for which it becomes significant. Even atoms have been considered as local 
systems |10[|11| . and a connection with the well known anomaly observed in 
the Pioneer satellites has been proposed [T^, although there is no foundation 
for attributing the Pioneer anomaly to the cosmic expansion [?] . Following the 
introduction of the dark energy concept to explain the present acceleration of 
the universe discovered with type la supernovae [13] and the realization that 
this dark energy could potentially take the form of phantom energy leading 
to a Big Rip singularity at a finite future [13], the effect of the cosmological 
expansion on local systems as the Big Rip is approached seems to go unques- 
tioned even though the persistence of cosmic effects on local systems in an 
adiabatic approximation in which the Hubble scale is much larger than the 
typical scales for local dynamics is often denied (see [5] for details). 

Two lines of approach to the problem of cosmological expansion versus local 
dynamics have been followed. The most common approach studies a Newto- 
nian gravitationally bound system, such as a binary stellar system embedded 
in an expanding FLRW universe and computes the effect of the cosmic ex- 
pansion as a perturbation of the local dynamics. The result for a particle of 
polar coordinates (r, and angular momentum L in the field of a mass M 
embedded in a FLRW universe with scale factor a{t) is given by the equations 
of motion [UlSJi]: 



The second approach, originally pursued with the McVittie solution [T] and 
the Einstein- Straus vacuole f^lSl used analytical solutions of the Einstein equa- 
tions to describe a local (spherically symmetric) inhomogeneity embedded in a 
FLRW universe (see [TB] for a review of inhomogeneous cosmological models) . 
The Lemaitre-Tolman-Bondi (LTB) class of solutions {e.g.^ [S]) describes a 
spherical inhomogeneity embedded in a dust-dominated (pressureless) FLRW 
universe. Here we propose, in the LTB class of solutions, a particularly clear 
example of a local system composed of two spherical shells with density merg- 
ing smoothly with the cosmic substratum, which is embedded in a FLRW 
universe. These shells have zero angular momentum and are kept at finite 
distance from each other by the cosmic expansion. Clearly, the latter has a 
non-negligible effect on the local dynamics of this system, independent of the 
strength of the local interaction. In fact, the effect of the cosmological expan- 
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sion on the local two-shell system persists even in the limit in which the ratio 
between local and cosmological densities is arbitrarily large. For simplicity, we 
consider a spatially flat FLRW background and we follow the notations of Ref. 
[T7] . The metric signature is — h ++ and units are used in which the speed of 
light in vacuo c and the gravitational constant G assume the value unity. 



2 An exact two body system in a cosmological background 

The LTB line element is HUZIIH] 



dt^ + [R'{t, r)f dr^ + (t, r) dH^ , 



(3) 



where r is a comoving radius, dH^ — diP -|-sin^ d dip^ is the metric on the unit 
2-sphere, 

2/3 

(4) 



is an areal radius, 



me{r) = Att dr' r po{r') , 



(5) 



is called the "Euclidean mass", and po{r) is the energy density on an initial 
hypersurface. The areal radius dH) is obtained as solution of the classical Bondi 
equation IHMIIIHIIS] 

me{r) 



(6) 



R{t,r) ■ 

A prime and a dot denote differentiation with respect to r and t, respectively. 

Let us consider, at the initial time, a density function that has two equal 
peaks located at ri and r2 (with r2 > ^i) 



Poir) ^ Pc + P* 
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where pc is a positive constant denoting the density of a cosmological back- 
ground far away from the density peaks, i.e., pa{r) w pc as r — ^ -l-oo, and 
is proportional to the peak density with p^jc? > pc. 

With the choice ([7]) for the initial density po{r), the Euclidean mass can 
be calculated explicitly as (see Figs. [U [5]) 
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Fig. 1 The initial density po(r) for the parameters choice pc = l,pt = 2, a? = 2,ri = 8, 
and r2 = 16 (in arbitrary units). 
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Fig. 2 The Euclidean mass m^ir). 
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This Euclidean mass is related to the Misner-Sharp-Hernandez mass and 
to the Hawking-Hayward quasi-local mass. The Misner-Sharp-Hernandez mass 
Mmsh is defined, for a spherically symmetric metric, by [TSlfT^ 



(9) 
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from which it fohows that 

r 

Mmsh = ^2nJ dspo{s)s' . (10) 



The first equality is independent of the particular form of the initial den- 
sity distribution po{r). Further, with spherical symmetry, the Misner-Sharp- 
Hernandez mass coincide^ with the Hawking-Hayward quasi-local mass [201 

With a slight abuse of terminology, the two density peaks at ri and r2 
will be called "shells". In actual fact, they are not sharp shells, but contin- 
uous distributions of matter more or less peaked according to the value of 
the parameter a appearing in (jT]), and with maximum density controlled by 
the parameter . This two- "shell" system mimicks the binary system of test 
particles chosen in many works to assess the effect of the cosmic expansion 
on local systems [1], but with two importance differences. First, there is zero 
angular momentum, no direction is preferred, and spherical symmetry holds 
because instead of a single test particle orbiting around a centre of force and 
subject to the cosmic expansion as a perturbation of this motion, we now have 
a spherical "shell" surrounding a smaller "shell" acting as a centre of force. 
Second, there is no test particle here: the dust with density p gravitates and 
curves spacetime. Instead of treating the cosmic expansion as a small pertur- 
bation of a weakly gravitating system, we have an exact solution describing 
at once the spacetime and the two-shells matter distribution, which merge 
smoothly with the cosmological background, and this is the virtue of our LTB 
example. No approximations are made and no adiabatic expansion in term of 
local and Hubble time scales is necessary. The element of proper radial dis- 
tance is (setting dt, dd, and dip equal to zero) dl = R'dr. Using the fact that 
dR = Rtdt + R'dr, we have dl = dR for dt = 0. Therefore, the proper radial 
distance between the two density peaks on a constant time slice of the LTB 
spacetime is 

h2{t) = i?2 - i?i = R{t, ra) - i?(i, ri) = 

This quantity can be regarded as the size of the two-shell system and corre- 
sponds, roughly speaking, to the radius of the binary system considered in 
previous works. However, I12 is now a proper distance, not a coordinate dis- 
tance (the use of misleading coordinate systems has led to coordinate-based 
statements and has marred many investigations of this problem) . At the time 
t = this proper radial distance coincides with the comoving coordinate dis- 
tance (r2 — ri), but it departs from it for t > 0. A plot of I12 is given in fig. [31 



1 See, e.g., Ref. [4]. 
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Fig. 3 The proper distance l\2(i) as a function of the comoving time. At t = 0, ii2 coincides 
with the comoving coordinate distance ■r2 — r'l . 



At large times i — !■ +oo, it is 
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i.e., the size I12 of the two-shell system is comoving and scales like the scale 
factor a{t) t^/^ of a dust-dominated FLRW universe, to which LTB space- 
times reduce for large r. This result matches the one for the Einstein-Straus 
vacuole [211511 22] in which a spherical central inhomogeneity is surrounded by 
a spherical vacuum, which is matched to a dust-dominated FLRW space. The 
proper radius of this matching surface is perfectly comoving [4] . 

We can now discuss the initial conditions in our model. The initial radial 
separation /i2(0) — R{0,r2) — i?(0,ri) = r2 — ri between the two "shells" can 
be taken arbitrarily. Following equation ([5]), their relative (radial) velocity at 
time t is 



Ut) -Ht.r.) ^Rit,n) = ±LteH - . (13) 



R{t,r2) y Rit,ri) 



At the initial time < = 0, it is 



The functions y^^mjpfjr and —y^rnjrjjr are plotted in fig. S) They are, 
respectively, always increasing and always decreasing, therefore the sign of the 
initial relative velocity ^12(0) of the two "shells" is conserved. Assuming that 
the relative initial velocity of the "shells" is not sufficiently large that the 
shells collide and cause shell singularities (then the formalism outlined above 
remains valid), the relative velocity of the "shells" I12 given by eq. (fT5)) tends 
to zero because R{t,r) oc t^/^. The "shells" then become comoving. 
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Fig. 4 The functions ^ 'me{r)/r and — ^ m,,{r)/r for the choice (0 of the initial density 
and the same parameter values of fig. ^ 

3 Discussion 

Generally speaking, the Hubble flow washes out the effects of a local over- 
density; this fact is well known in FLRW universes and it occurs also in LTB 
spacetimes. The following argument is a bit naive but nevertheless carries some 
weight: consider a function /(r) and its radial gradient ^ with respect to the 
physical (areal) radius R (the only spatial gradient allowed in the presence of 
spherical symmetry). It is straightforward to see that 



and, therefore, df /dR — > as i — +00. The evolution of our LTB model is 
in agreement with this general feature. 

In our LTB model, in which the density contrast is smoothed out by the 
cosmic expansion, the energy density on an initial hypersurface can be chosen 
at will and the expression ([T]) describes two local concentrations of matter 
("shells") merging smoothly with the cosmological background. These two 
"shells" play a role analogous to that of a binary system composed of two 
test particles embedded in a FLRW background universe so often considered 
in studies of the effects of the cosmological expansion on local systems. We 
have chosen the peaks of the initial density distribution to be of equal height 
for simplicity but this assumption is not necessary: the ratio of their heights 
can be chosen at will, in the same way that the masses of a binary system 
subject to the effect of the cosmological expansion can be chosen in any ra- 
tio. The LTB example that we presented has two clear advantages over the 
binary: first, it is a fully gravitating configuration, not a system of test par- 
ticles and one does not need to resort to approximations in order to embed 
this system into a cosmological background. This is an exact solution of the 
Einstein equations. Second, due to the non linearity of general relativity, when 
one wants to consider cosmic expansion versus local dynamics, it is a ■priori 
difficult or impossible to split a solution into a background plus an inhomo- 
geneity (apart from perturbative regimes, which are nevertheless subject to the 
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notorious gauge-dependence problems and require the use of gauge-invariant 
formalisms): this is exactly what the LTB metric (or the McVittie metric [U 
[nm]) does for us. 

It is reasonable to regard the proper radial distance I12 = R2 — R1 between 
the two "shells" as the physical size of the local system. This quantity depends 
on time and it is straightforward to see that this distance is comoving with 
the cosmic dust of the FLRW background. In the LTB example proposed, 
the local system is affected by the cosmic expansion, regardless of its spatial 
size. One can take the limit — > in which the system becomes a test 
fluid, or the opposite limit p*/a^ » pc in which the local system is strongly 
gravitating. The initial distance R2 ^ Ri = r2 — ri between the two density 
peaks can be adjusted at will in comparison with the Hubble radius H^^ of 
the FLRW background. Any way these parameters are varied, the result is 
always the same: the local system is stretched by the cosmological expansion 
and participates in it. This fact shows that, assuming the FLRW metric to 
describe the universe down to small scales, it is not true that systems below 
a certain spatial scale are unaffected while only larger ones partake into the 
expansion, as suggested by many authors (see [5l|4] for references) . By choosing 
the LTB metric as an example, we were bound to find this result; indeed, it is 
built into the LTB metric itself. One may, therefore, question the validity of our 
argument. The point is that a solution with these properties does exist, it is a 
perfectly legitimate example, and one of the rare exact solutions of the Einstein 
equations in which a clear answer to the puzzle of cosmic expansion versus 
local dynamics can be obtained. The answer matches the results obtained with 
the Einstein- Straus vacuole McVittie generahzed McVittie 

solutions [5], and LTB black holes [23]. From the physical point of view, a 
relevant issue is whether the two-shell system is stable, but this question can 
also be posed for any LTB metric. 

One potential problem with the LTB model proposed here consists of the 
possible presence of shell-crossing singularities in the t < region of this 
spacetime. Although these singularities are not mandatory in LTB models 
(conditions for their avoidance are given in Ref. [28]), they are to be expected 
in general. The prevailing opinion about shell-crossing singularities is that they 
are an artifact of the dust equation of state chosen and should disappear when, 
more realistically, some pressure is introduced (see, e.g., Ref. [IS]) although, to 
be honest, we do not know of a rigorous mathematical proof of this statement 
to date. Overall, we regard the possibility of shell-crossing singularities as a 
non-essential feature of LTB models and we choose to focus on the region of 
LTB spacetime which is free of these singularities. 

Finally, we could make the two shells of our model infinitely thin and, 
in this case, they would have to be matched to their surroundings by using 
the Darmois-Israel junction conditions [25| . However, this construction is not 
necessary and by continuity of the property presented here, one expects that a 
system composed of two zero-thickness shells will also be comoving. The case 
of a single shell matched through the junction conditions to an expanding 
FLRW universe, and hosting a wormhole, was considered in Refs. |^[?7] with 
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the result that, even if the shell initially has a non- vanishing velocity relative 
to the cosmic substratum, it eventually becomes comoving with it. The answer 
to the question of cosmic expansion versus local dynamics seems to be that 
expansion wins and local systems go with the (Hubble) flow. 
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